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Abstract. The idea of the meson cloud approach in the chiral quark model has been 
extended to include gluon cloud in order to achieve the parton densities in the nuclcon, 
based on the constitute quark framework. The splitting function of the quark to the 
quark-meson and quark-gluon at low Q 2 value are used to obtain parton densities in 
the constituent quark. The phenomcnological constituent model is employed to extract 
the parton distributions in the proton at low Q 2 value. Since we have access to the 
parton densities at low Q 2 , we are able to obtain F2(x,Q 2 ) structure function at low 
Q 2 value. The result is in good agreement with available experimental data and some 
theoretical models. To confirm the validity of our calculations, the fraction of total 
momentum of proton which is carried by gluon at high Q 2 and also the Gottfried sum 
rule are computed. The results are in good agreement with what are expected. 



1. Introduction 

In hard scattering events, mesons and baryons can be viewed as bound states built up 
from partonic constituents, i.e. quarks and gluons. This picture changes at low energies, 
where hadronic effects play a more prominent role in the non-perturbative structure of 
hadrons. One particular importance is the pion cloud effects which e.g. have a direct 
impact on the spin structure of the proton pp. Thus they need to be incorporated in 
bound-state calculations aiming at a realistic description of mesons and baryons. Pion 
effects on the quark propagation are important for several reasons. They account for (at 
least part of the) pion cloud effects in baryons and mesons. Furthermore they allow for 
the possibility of hadronic intermediate states in bound state calculations and therefore 
generate the finite width of meson spectral functions. 
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Dynamical chiral symmetry breaking is one of the most important properties of low 
energy QCD. The breaking pattern has profound impact for phenomenological quanti- 
ties, as e.g. the appearance of the pseudoscalar Goldstone bosons in the chiral limit of 
QCD and the non-degeneracy of chiral partners. Chiral perturbation theory [2J [3] de- 
scribes these effects very efficiently on the level of hadrons but has nothing to say about 
the underlying structure of the full theory. The interplay between the fundamental quark 
and gluon degrees of freedom and the resulting bound states are also particularly inter- 
esting. In full QCD there are hadronic contributions to the fully dressed quark-gluon 
interaction. These effects are generated by the inclusion of dynamical sea quarks in the 
quark-gluon interaction and are therefore only present in unquenched QCD. The quark- 
gluon vertex is also an important ingredient into the quark-antiquark interaction that is 
responsible for the formation and properties of bound states. On a perturbative level, 
the quark-gluon vertex has been studied in detail in arbitrary gauge and dimensions in 
[I] . However nonperturbative properties of this vertex are also still under tense scrutiny. 

At low energies, the idea that baryons are made up of three constituent quarks 
and mesons of a (constituent) quark-antiquark pair , the naive quark model scenario, 
accounts for a large number of experimental facts . The quest for a relation between the 
two regimes, i.e. between the current quarks of the theory and the constituent quarks 
of the model has an old history and, in recent years, this search has been the subject of 
a considerable research effort . The fundamental problem one would like to understand 
is how confinement, i.e. the apparent absence of color charges and dynamics in hadron 
physics, is realized. Detailed quark models of hadron structure based on the constituent 
quark concept have been defined in order to explain low energy properties To pro- 
ceed from these models to the asymptotic regime, where deep inelastic scattering (DIS) 
takes place, a hadronic scale is associated to the model calculations. The experimental 
conditions are reached by projecting the leading twist component of the observable and 
evolving according to perturbative QCD. The procedure describes successfully the gross 
features of the DIS results. It was long ago, at the time that QCD was being proposed, 
that a procedure, hereafter called ACMP(Altarelli, Cabibo, Miani, Petronzio) [6], was 
developed to understand the relation between the constituent quarks and the partons. 
In this approach, constituent quarks are complex objects, made up of point-like partons 
(current quarks, antiquarks and gluons), interacting by a residual interaction described 
by a quark model. The hadron structure functions are obtained as a convolution of the 
constituent quark wave function with the constituent quark structure function. This 
procedure has been recently reviewed to estimate the structure function of the pion 
with success. In the ACMP approach, each constituent quark is dressed by a neutral 
cloud of quark-antiquark pairs and gluons, thus, this scenario supports the confinement 
mechanism. A few years earlier a second approach had been developed [7] , in which the 
proton is assumed to be made out of three valence quarks plus a neutral core of quark- 
antiquark pairs and gluons, very much in the spirit of recent developments along the 
Manohar-Georgi model [5]. This duality of approaches has to do, in modern language, 
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with the implementation of Chiral Symmetry Breaking(CSB). The naive models do not 
contain spontaneous CSB and this phenomenon has to be implemented if they are to 
represent QCD at low energies. 

The effective chiral quark theory [8] may provide an alternative explanation to that 
of the traditional meson cloud approach [I9j. In this theory, the relevant degrees of 
freedom are constituent quarks, gluons and Goldstone bosons. The chiral quark model 
(xQM) includes both gluon and pion exchange between constituent quarks together 
with corresponding exchange currents. The relevant degrees of freedom and the related 
question whether the pions couple effectively to the nucleon or to the constituent quarks 
is extensively discussed in Refs. [9], [10] . It is necessary to study the consequences of these 
different scenarios in a broad range of physical processes to assess their validity. 

On the other hand, the xQM can be used to study the flavor structure of the 
constituent quark model and the nucleon within the conventional mesonic cloud pic- 
ture. Using this model the effects of SU(3)f symmetry breaking can be discussed [TT] . 
The implications of the Gottfried sum rule (GSR) violation for the A-n mass splitting 
were also considered in Ref.JTT]. At low energy resolution scale the constituent quark 
picture successfully describes hadronic structure functions. The sea quark and gluonic 
degrees of freedom are assumed to be absorbed into constituent quarks to be considered 
as quasi-particles [12] . A relation between the two regimes of hadron structure function 
description; i.e. the chiral quark and the constituent quark models, has a considerable 
significance which has been investigated widely in the literature, and has attracted much 
attention in recent years [T3] . 

As should be noticed, the main ingredients of this paper are two subjects. In con- 
tinuation of our previous work [2] we add the gluon cloud to the xQM while we use an 
effective lagrangian at low Q 2 values. We resort to a constituent quark model to extract 
parton densities inside the proton. Since the gluon densities are also at our access, we 
are able to calculate F 2 structure function for the proton at the NLO approximation. 

The organization of the paper is as follows: In section 2 we introduce Quark-meson 
and in similar way quark-gluon vertex function based on nucleonic Sullivan deep- inelastic 
scattering [15] . Section 3 is allocated to xQM and the constituent quark distribution is 
obtained. In addition we consider the gluon cloud in the constituent quark. Therefore 
we can calculate the gluon distribution function inside the proton. This contribution has 
not been considered in Ref.[TT]. In this section we also discuss about phenomenological 
valon model [16] to extract valence distribution in the meson. These distributions are 
required to obtain mesonic anti-quark contribution in the constituent quarks. In section 
4 we present our result for F 2 structure function at the NLO approximation, using the 
parton densities in the proton which are extracted from xQM. To confirm the validity 
of our calculation, we evolve the gluon distribution to high Q 2 to get the momentum 



LowQ 2 proton structure function, using gluon and pseudoscalar meson clouds in the constituent quark frame 



fraction of proton which is carried by gluon. Using the antisymmetric property of sea 
quark densities which is resulted from xQM, we also calculate the Gottfired sum rule 
to test again the validity of our calculations. The conclusion is given in section 5. 

2. Chiral quark model and effective quark-meson and quark-gluon 
interactions 

This model was introduced by Georgi and Manohar [8] in order to incorporate the chiral 
symmetry of QCD into the successful features of the constituent quark model [T7]. The 
prime assumption of the model is the possible realization of an effective Lagrangian 
between the scale of chiral symmetry breaking A x and confinement scale Aqcd- The 
dynamical degrees of freedom here are constituent quarks, pseudoscalar mesons and 
arguably gluons. The respecting Lagrangian can be written as 

Ceff = iQl^ + ig s G»)Q - ^Qd,UY^ 5 Q 

Jit 

-M q QQ + ^traciD^UD^) - \g^G^ . (1) 

The matrix U contains the pseudoscalar mesons and Q stands for the constituent quark. 
G^ v is the strength field tensor defined as following 

G""'" = d li G u,a - d u G^ a - if abc G^G v c (2) 

where G^ is the gluon field and f 's are the structure constants. The covariant derivative, 
D^, is associated to the chiral symmetry of QCD in flavour space. 

Mq, f n and Qa are the constituent quark mass, the pion decay constant and the 
axial- vector constant, respectively. The strong running coupling, g s , has to be consid- 
ered in some energy below A x and we then take it as a constant. 



2.1. Quark-meson effective vertex function 

The effective chiral quark model [8] is applied in order to study the pseudoscalar meson 
clouds in the constituent quarks [UJ. In Ref.pXj it is found that the pionic clouds can 
explain the violation of the Gottfried sum rule and introduce an enhancement on the 
non-perturbative effects of the sea quark pairs. To this end the pion-quark splitting 
function is introduced in analogy to the nucleonic Sullivan deep-inelastic scattering and 
expressed in [15] as: 

r f i 2 \ 9q->mq' 1 I n i j 2 \ 1 2 

Jq^mq'{Xm,k±) — — 2 7z r \ < ^q^mq>{Xm,k±)\ 

\m l xm(1 - x M ) 

((1 - x M )m Q - m Ql f + k\ 

(l-x M )(ml-M MQ ,y ' [6) 

where % is the longitudinal (light cone) momentum fraction of the constituent quark 
for the Goldstone boson and k± is the perpendicular momentum of the quark Q' . The 
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g is the effective coupling constant of pion- constituent quark: 

2 _ 9 a (™Q + m Q >) 2 
9q^mq> p 4 ' V 1 ) 

where is the axial vector coupling constant of the constituent quark which is equal 
to one. We take mi = vtlqi = ^f- = 313 MeV for the light up and down quarks and 
m s = rriQt = ms — m^ + mi = 567 MeV for the strange quarks. 

2 

MQ' 

2 m 2 M + k 2 ± m QI + k\ 

M mq> = + -j — • (5) 

x M 1 — X M 



M\,qi is the invariant mass squared of MQ' system which is defined as: 



The Gq^mq 1 is a vertex function or phenomenological form factor for which we adopt 
the exponential form: 

f m 2 Q - M 2 MQ ,(x M ,kl) \ 
Gq^mq' = exp ^ ^ j . (6) 

A x is a cut-off parameter which can be taken equal for all fluctuations involving 
pseudoscalar or vector mesons. The integration of the splitting function over k± 

[fQ->MQ'(% M ) = jo°° fQ^MQ>(x M ,kr) <^t] and tnen over x m and finally summing over 
the intermediate quarks (Q') yields: 



Pm/q 



\ a M/q\ 2 = Y] / f Q ^ M Qi(x M )dx M , (7) 

r^l JO 



which is the probability of finding a Goldstone boson M in the constituent quark Q. 



2.2. Quark-gluon effective vertex function 

Fairly gluon distributions can be obtained by dressing quarks with gluons in the 
nonperturbative regime with massive effective gluons {m e J* ) and frozen running a s . 
Rather heavy effective gluons m e g " > 0.4 GeV and small a s < 0.5 are required in order 
to limit the momentum carried by quarks to approximately what is required by the 
phenomenology [18]. Now in order to include the gluon clouds in the constituent quarks 
we need to put the line of analogy some further, in the sense that an almost the same 
form of splitting function is regarded for the gluon-quark interaction as that for the 
quark-meson interaction. The main differences stand on two parts. The first one is the 
quark- meson coupling constant, that we replace it with the strong coupling constant at 
some low energy. Secondly, we need to know the relevant vertex function for the quark- 
gluon interaction. The vertex function encodes the extended structure of the gluon 
and the constituent quarks. The extraction of the vertex function is rather difficult 
since it incorporates the non-perturbative effects. However, in a series of recent studies 
[20J [21] , t ne authors have calculated the non-perturbative corrections to the quark-gluon 
vertex in the framework of the Dyson- Schwinger and Bethe-Salpeter equation. Their 
predictions for the light meson properties seems satisfactory [20]. On the other hand, we 
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have found out that our ansatz for the quark-gluon vertex which is assumed similar to 
the quark-meson vertex, has qualitatively the same momentum behavior. 

Consequently, we have the quark-gluon fluctuations which tends to the following 
splitting function: 

fQ^9Q'(Xg,k 2 ± ) = ^ - - ^ - - ^ \GQ^gQ>{Xg,k]_)\ 

((1 - %g) m Q ~ m Q'Y + fcj f M 
{l-x g ){ml-Ml Ql f ' 1 ] 

where x g is the longitudinal (light cone) momentum fraction of the constituent quark 
for the gluon and k± is the perpendicular momentum of the quark Q' . 

The integration of the quark-gluon splitting function over k± and then over x g and 
finally summing over the intermediate quarks (Q f ) yields: 

P g/Q = K/Ql 2 = J2 I fQ^gQ'( x g ) dx g , ( 9 ) 



Q' 



o 



3. Constituent quark distribution function in the chiral quark model 

The constituent quark Fock-state \Q), can be expressed in terms of a series of light-cone 
Fock-states: 

\Q) = Vz\q) + J2*B/ Q \q',B), (10) 

q' 

where \q) is the "bare" but massive state, yZ denotes the renormalization factor for 
a "bare" constituent quark and |ob/q| 2 are probabilities to find Goldstone bosons and 
gluon distribution in the constituent quark states. Then the dressed u- and d-quark 
Fock-states are: 

\U) = \[Z\u) + a n0/u \u, vr°) + a n+/u \d, vr + ) + a k+/u \s, K + ) 

+ a g /u\u, g) + ■ ■ ■ , (11) 

\D) = \/Z\d, ) + ^ a n o /D \d, 7T°) + ^1 a„-/ D \u, vr") + a k o /D \s, K°) 

+ a g/D \d,g) H . (12) 

The above expressions can be depicted as in the following graphs: 

7T° 7T+ K+ 9 



U u uuuud u usuuuu 



7T° 7^ 



K° 9 



D d dddd u dd s dddd 
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In these graphs the thick lines indicate the quark propagators, dashed lines are 
meson fields and wiggly curves stand for the gluons, respectively. 



We consider the nucleon to be a bound state of three constituent quarks (U and 
D). The quark distributions in the constituent quark, at some QCD initial scale, can 
be written as: 



u n (x) 


= uP(x) + % \x) 


+ u^\x) + u^(x), 


(13) 


u D (x) 






(14) 


d D (x) 


= d ( °\x) + d ( $(x) 


+ d ( £\x) + d%\x), 


(15) 


du(x) 


= d { £\x) + d { u\x) 


1 


(16) 


SU{D){X) 


= S U(D) + S U(D)( X . 


), 


(17) 


g v (x) 


= u i §\l-x) 




(18) 


9d(x) 


= dg>(l-s). 




(19) 



%( 


» 


-W 

= % 


= d D { 


x ) 


— a D . 


U D 


(x) 


-W 
= u K D > 


= d v ( 


x) 


— Ujj , 






_ J K ) 
— b u 


= s D [ 


» 


-s (K) 






~(k) 
= i% 


- d {k) 

— u D 







The anti-quark distributions become: 

(20) 
(21) 
(22) 
(23) 

It should be noted that for sea quark densities we have = q~Q ' . The superscripts 
denoted with (0) correspond to the "bare" quark distributions, and those denoted with 
(i and g) to the intermediate quark distributions associated with mesons and gluons 
respectively, and those denoted with (it) originate from mesons (pions). 
The bare quark distribution in the constituent quarks has the form: 

(o)/^ _ JP), 



U\j'[X 



d%\x)=\l-J^P c/Q \5(x-l), (24) 

where these distributions play the role of the valance quark distributions inside the 
constituent quarks. In Eq. ( )24|) . Pc/Q refers to the probability of finding a Goldstone 
boson and gluon in the constituent quark Q. So back to Eq.liTUl). we have Pc/Q = \cib/q\ 2 - 

The intermediate quark distribution function in the constituent quark is calculated 
from the meson splitting function: 

4\x) = d { ${x) = \u /Q {l-x), (25) 
«g ) (x) = d«(x) = 1/^(1 -x), (26) 



x)=f K/Q (l-x), (27) 
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where / m is the total splitting function of the constituent quark and is defined as 

Q 

% = Yl fQ->MQ' ■ ( 28 ) 
Q' 

Mesonic anti-quark contributions in the constituent quarks which finally yield us sea 
quark distributions inside the constituent quarks, are given by the equations: 

v$\x) = ±I K {x), u<S>(x) = ^I K {x), s\j K \x) = I K (x), (29) 

where 

Im{x) = J SmJq{v)Qm y • (30) 

Here qu{y) denotes the valence quark distribution of the meson. These valence 
distributions are required to extract sea quark densities in the constituent quarks of the 
proton. They will be obtained in the next section, using the phenomenological valon 
model. 

Using the constituent quark model for the proton , the parton densities at low Q 2 
value, for instance 0.5 GeV 2 can be obtained. More details regarding the employed 
constituent model shall be explained in section. 4. 

3.1. Valon model 

According to the valon model [16J, a valon is a dressed valence quark so that there is a 

one-to-one identification of a valon with the associated valence quark as probed at high 

Q 2 . In this model a meson, for instance, is a bound state of two valons. They contribute 

independently in an inclusive hard collision with a Q 2 -dependence that can be calculated 

in QCD at high Q 2 . The valon picture suggests that the structure function of a meson 

involves a convolution of two distributions: the valon distribution in the meson and the 

structure function for each valon, so that one has 

i 



F 2 M (x,Q 2 ) = / dyG v/M (y)F2(z = -,Q 2 ) , (31) 

v Jx y 

where the summation is over the two valons. Here F 2 (z, Q 2 ) is the meson structure 
function, F 2 is the corresponding structure function of a v valon, and G v /m{v) indicates 
the probability for the v valon to have momentum fraction y in the meson. We shall 
assume that the two valons carry all the momentum of the meson. 
We assume the following simple form for the exclusive valon distribution inside the 
mesons which facilitates the phenomenological analysis, 

G v ( yi ,y 2 ) = g(yi) p y q 2 5(yi + y2 - 1), (32) 
where p and q are two free parameters and is the momentum fraction of the i'th valon. 
The U and D type inclusive valon distributions can be obtained by integration over the 
specified variable, 

G Vl (y)= f dy 2 G v (y,y 2 )=gy p (l- y y, (33) 
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G V2 {y)= f dy 1 G v (y 1 ,y) = gy ll (l-yy. (34) 



The normalization parameter g has been fixed by requiring 

G Vl (y)dy= G V2 (y)dy = 1, (35) 
o Jo 

and is given by g = B ^ p+ \ , where B(m,n) is the Euler-beta function. 
Consequently, we will get the following inclusive valon distributions for mesons: 

The dirac delta function, 5(yi +y 2 — l), automatically ensures the momentum sum rule: 

yG Vl dy+ yG V2 dy = 1 , (37) 
'o Jo 

Mellin transformation from Eq. (I3T)|) . will yield the moment distributions for valons 
[13]. The moments of quark and gluon distributions at any energy scale inside the 
meson are obtained by multiplying the valon moments with the appropriate moments 
of singlet, non-singlet and gluon sectors. Using the inverse Mellin transformation in the 
parameterized form as described in [22] and fitting over the available experimental data, 
the valence quark densities inside the mesons will be obtained. 

4. Results and discussions 

Using the xQM we are able to extract the valence, sea and gluon densities inside 
the constituent quarks. To access to parton densities inside the proton, we employ 
a constituent model. We need the quark distribution in a proton, gjv(x), which 
can be obtained using the convolution of the corresponding quark distributions in 
the constituent quark (qu,D(x/y)) with the light-cone momentum distribution of the 
constituent quark in the proton (U(y),D(y)), so that : 



q N {x) 



2U(y)q u (-)+D(y)q D (- 



dy 



(38) 

y 



mj \y 

Eq. (J38l) is the basis for the constituent quark model in which we can use to obtain the 
quark densities in a proton. In our calculations, U(y) and D(y) are parameterized as: 

U(y) = 0(a Q + l,f3 Q + 1)^(1 - y) B ® 

D(y) = 0{l Q + hv Q + 1)2/^(1 - y) VQ (39) 
where B(ag + 1,{3q + 1) and Bi^q + 1,T]q + 1) are the Euler beta functions. 
The normalization coefficient B ^ a + { Bq+1 ^ and B ^ Q+ [ VQ +i) anc ^ fi na Uy the unknown 



parameters which exist in Eq. ( |39|) have been fixed by requiring the number sum rule 
for valance quark densities inside the proton and also the momentum sum rule for the 
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X 



Figure 1. Sea quark distribution in the proton at low Q 2 = 0.5 GeV 2 , based on xQM ■ 
The anti-symmetric of the sea quark distributions are obvious. 

parton densities inside the proton. 

Requiring sum rules, we will get the following numerical values for the parameters 
associated with the constituent quark distributions : 



U(y) 


a v 


Pu 




-0.40 


0.57 


D(y) 


lu 


Vu 




-0.39 


0.55 



The obtained densities are at low Q 2 =0.5 GeV 2 which is in correspond to the chosen 
value for A^l.26 GeV in our calculations according to the model A of Ref . |H| . 
The results for sea and gluon densities in the proton are depicted in Fig.l and 
Fig. 2 respectively. The asymmetry of sea quark densities are obvious as we expected 
from xQM. By using the parton densities, the F 2 structure function at the NLO 
approximation which is defined by 

F* p (x, Q 2 ) = xJ2 Q 2 ) + fa, Q 2 ) + x KM + q) + 2 C£ >2 g}} (40) 

q 

can be calculated. In Eq.f l40p q, q and g refer respectively to quark and gluon distri- 
bution inside the proton and C terms are Wilson coefficients which are defined in [18]. 
In Fig.3 the F 2 structure function for the proton at low Q 2 = 0.5 GeV 2 is shown and 
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Figure 2. Gluon distribution in the proton at low Q 2 = 0.5 GeV 2 , based on xQM ■ A 
comparison with GRSV has also been done. 

compared with the GRSV model [18J. A comparison with available experimental data 
[19] has also been done. The agreement is well. Since we have access to gluon distri- 
bution, to confirm the validity of calculation at low Q 2 value, we can evolve it to high 
Q 2 and calculate the fraction momentum of proton carried by gluon. In this regard 
we got 41.2% which is what we expect. Authors in [T8] claimed that their extracted 
gluon distribution carries about 50% of total momentum of the proton. This is more 
than what we got. On viewing at Fig. 2 which indicates gluon distribution in the proton 
at low Q 2 , it is predictable that we reach to lower presentation of fraction momentum 
at evolved Q 2 value in our model. A justifiable reason for this difference backs to this 
reality that the results in [18] is based on a global fit while we employ xQM to obtain 
the gluon contribution in our calculation. 

We use the following relation [TTJ 

f 1 dr 1 2 f l - 14 

Sg = I - F?(x)) — = 3 + 3 ^ " d{x))dx =3-5^, (41) 

to obtain the Gottfried sum rule (GSR). The numerical value which is obtained by this 
model is 0.2339 which is very near to quoted experimental value 0.235 ± 0.026 by NMC 
group [231 121] ■ Once again the validity of the calculation using xQM at low Q 2 value 
is confirmed. 
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Figure 3. Analytical result for F% structure function at LO and NLO approximation. 
The result have been compared with GRSV model [TSJ and available experimental data 



5. Conclusion 

The flavor structure of the nucleon in the effective chiral quark model (xQM) has 
been studied. In this model the Goldstone bosons couple directly to the constituent 
quarks. This idea has been extended to include gluon cloud in the xQM at low Q 2 
value. Consequently we could obtain an improvement result for sea quark density and 
also calculated the gluon distribution directly in xQM while it has not been done in 
previous works [TTJ [H] . To obtain the sea quark densities inside the constituent quark, 
we needed valence quark distributions of the meson . We got these valence distributions 
using the phenomenological valon model [16J. Furthermore the sea quark densities in 
the proton have been obtained, by convoluting the required distributions in the used 
constituent quark model. This yielded a result in which the sea quark densities in the 
proton are un-symmetrized. By accessing to the other partons in the proton, we could 
calculate the F% structure function at Q 2 = 0.5 GeV 2 which confirms the anticipated 
result of the model. For more validity of the model the fraction momentum of proton 
which is carried by gluon at Q 2 = 15 GeV 2 has been calculated. The numerical result 
which was obtained for this fraction at Q 2 = 15 GeV 2 and also the numerical value for 
GSR are very close to what are expected. However the obtained results in this paper 
are satisfactory but one can use different vertex function, quoted in [25J . By comparing 
the results, one can choose the best candidate to consider meson and gluon clouds at 
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low values Q 2 . This can be done as a new research job in future. Further suggestion 
is to consider valence quark density inside the constituent quark with a gaussian form 
rather than the Dirac Delta as in Eq. In this case we would expect to achieve the 

number sum rules for the both constituent and valence quarks in a more straightforward 
way. This will also be a scientific challenge for a research task in future. 
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